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\ We study signals for the production of superparticles at the Tevatron in supergravity 

' scenarios based on the Grand Unified group SO (10). The breaking of this group 

introduces extra contributions to the masses of all scalars, described by a single new 
Qh! parameter. We find that varying this parameter can considerably change the size 

of various expected signals studied in the literature, with different numbers of jets 
and/or charged leptons in the final state. The ratios of these signal can thus serve as a 
diagnostic to detect or constrain deviations from the much-studied scenario where all 
scalar masses are universal at the GUT scale. Moreover, under favorable circumstances 
^ ■ some of these signals, and/or new signals involving hard b— jets, should be observable 

at the next run of the Tevatron collider even if the average scalar mass lies well above 
the gluino mass. 
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1 Introduction 



The search for Supersymmetry (SUSY) [1], the most attractive candidate for physics beyond 
the Standard Model (SM), has top priority in present day high energy experiments. From 
the non-observation of sparticles at various colliders only lower limits on their masses have 
been obtained so far. 

From the model building point of view the limit M~± > 91 GeV on the mass of the 
lighter chargino xf from LEP2 is of special importance. This limit, in conjunction with the 
popular assumption of gaugino mass unification |lj with a common gaugino mass m^j at 
the GUT scale (Mq), implies that the gluino (g) must be rather heavy with m g ^ 300 GeV. 

The limit on M~± becomes much weaker if the xf happens to be approximately de- 
generate with either the lightest supersymmetric particle (LSP) assumed to be the lightest 
neutralino (x?), or with a sneutrino (u). The first possibility is strongly disfavored in models 
with gaugino mass unification and radiative electroweak symmetry breaking (see below). 
The second possibility occurs, e.g., in mSUGRA models JIJ for small values of the scalar soft 
breaking mass m . However, a recent result from the D0 collaboration || ruled out m g < 
300 GeV for small m from direct squark and gluino searches in the multi-jet plus missing 
Et channel. 

Like most other limits on squark and gluino production at hadron colliders, this bound 
holds if the squarks are degenerate with each other at Mq with a common mass mo- They are 
then also approximately degenerate at the weak scale with average mass m 9 -. Under the same 
assumption, for much heavier squarks (m, 3> nig) only m g % 200 GeV can be excluded from 
present Tevatron data. In this case, even gluino searches after the main injector upgrade 
of the Tevatron (referred to hereafter as Tevatron Run II), with an estimated integrated 
luminosity of 2 fb -1 at \/s = 2 TeV, are not likely to improve the indirect limit on the gluino 
mass from LEP2 searches, if gaugino mass unification holds. A further luminosity upgrade, 
amounting to an integrated luminosity ~ 25 fb _1 at y/s = 2 TeV (hereafter referred to as 
TeV33), could reach somewhat higher, but not dramatically different, masses through direct 
searches for gluinos and squarks |4jj. 

As already mentioned, these estimated search limits have been derived under certain 
assumptions, like gaugino mass unification and degeneracy of squarks at Mq- Such assump- 
tions about physics at high scales, although attractive, may ultimately prove wrong. It is 
therefore important to reexamine the search prospects at the Tevatron and its various up- 
grades if some or all of these assumptions are relaxed. If accompanied by suitable theoretical 
guidelines, this can be accomplished while still avoiding an unmanageable proliferation in 
the number of unknown parameters. Moreover, such an analysis indicates the feasibility of 
obtaining glimpses of the SUSY breaking pattern at a high scale using the Tevatron data. 
Finally, and most importantly, such analyses are needed to assess the robustness of the SUSY 
search strategies currently employed by Tevatron experiments. 

One possible avenue would be to relax the gaugino mass unification condition. This 
has been attempted by several groups ||. However, at least in the framework of SUSY 
GUT's, the assumption of gaugino mass unification appears quite natural, since it simply 
follows if the GUT symmetry is respected by the SUSY breaking mechanism at a high scale, 
irrespective of the specific choice of the GUT group. In this paper we, therefore, assume 
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that gaugino masses do indeed unify. 

The assumption that all soft breaking masses for the scalars have a universal value m at 
scale Mq (referred to hereafter as the conventional scenario), which in turn predicts approx- 
imately degenerate squarks at the weak scale, is more model dependent. In the framework 
of the popular N — 1 SUGRA GUT models 0, where SUSY breaking in the hidden sec- 
tor is transmitted to the observable sector by gravitational-strength interactions alone, two 
assumptions are required to obtain a degenerate scalar spectrum at scale Mq. First, the 
Kahler metric has to be the same for all chiral superfields. This ansatz drastically reduces 
the large number of parameters that may otherwise be present in the Minimal Supersym- 
metric Standard Model (MSSM), making the model much more predictive. It also avoids 
problems with flavor changing neutral currents caused by sparticle loops. This assumption 
guarantees degenerate sfermions at the scale A at which SUSY breaking is transmitted to 
the visible sector. By taking sfermions to be degenerate at scale Mq one therefore implicitly 
makes the second assumption that A ~ Mq ~ 2 • 10 16 GeV. However, from the supergravity 
point of view, it is more natural to choose A to be near the Planck scale Mp ~ 2.4 • 10 18 
GeV. Non-universality in scalar masses at scale Mq can then arise due to several reasons: 
i) The running of scalar masses between Mp and Mq |7j may lead to non-universality at Mq. 
Within the framework of an SO (10) SUSY GUT, however, the resulting non-universality is 
expected to be negligible for the first two generations of squarks and sleptons with small 
Yukawa couplings. In contrast, the soft breaking masses of the SU{2) doublet Higgs bosons 
responsible for electroweak symmetry breaking can be changed significantly, since they have 
to couple to super-heavy GUT fields in order to generate very large masses for their GUT 
partners (doublet-triplet splitting); moreover, these Higgs bosons are usually assumed to 
reside in a different representation of SO (10) than the sfermions do, and hence are renor- 
malized differently by gaugino loop diagrams. The situation for third generation sfermions 
is intermediate between these two extremes: they reside in the same representation as first 
and second generation sfermions do, but at least some of them have significant Yukawa cou- 
plings. Thus only the first two generations of sleptons and squarks can safely be considered 
to have a common mass at Mq, which we define to be m . The soft breaking masses of third 
generation sfermions are likely to be somewhat, but not very much, smaller, while those 
for the Higgs bosons can be either larger or smaller than m . Fortunately, the physics of 
the Higgs sector is relatively unimportant for direct squark and gluino searches at Tevatron 
energies. In our study we will neglect the (probably not very large) effect of the running 
between A and Mq on third generation soft breaking terms. 

ii) It has been known for some time that if the rank of a GUT group (or some symmetry 
group broken at an intermediate scale) is reduced by spontaneous symmetry breaking, one 
may obtain D— term contributions to scalar masses || [J. These will in general differ for 
different members of the same GUT multiplet, leading to non- universal squark and slepton 
masses at the symmetry breaking scale. The size of these new contributions can be com- 
parable to rriQ. Note that these nonuniversal terms are generation-independent, so that no 
additional problems due to flavor changing neutral currents arise. 

It has already been emphasized in |H| that right-handed down-type squarks (dp, Sp, bp, 
which are degenerate in mass) may be considerably lighter than other species of squarks 
(uL,dp or up) due to mechanism (ii). A specific model in which this may happen is an 
SO (10) SUSY GUT, breaking down to the SM either directly or via the Pati-Salam group 
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fTT| at an intermediate scale 0. In addition, the L-type (left-handed) or R-type (right- 
handed) sleptons may also turn out to be considerably lighter, which may yield many novel 
experimental signatures [[II], jl2f . The phenomenological consequences of D— terms at high 
scales have also been studied in the first paper of ||. 

In the present analysis we shall work within the framework of a SO( 10) SUSY GUT 
breaking down to the SM directly. In this model nonuniversality in the squark and slepton 
masses at a high scale can be parametrized by only one extra parameter D (compared to the 
conventional scenario) 0. 50(10) has several attractive features: it allows to fit a complete 
generation of (s)quarks and (s)leptons into a single representation of the group. Moreover, 
it includes right-handed neutrinos, and thus naturally allows to implement the see-saw 
mechanism |13| to explain the small but nonvanishing neutrino masses that seem to be needed 
to explain recent data on neutrino oscillations. Direct breaking to SU(3) x SU(2) x U(l) 
is perhaps more natural, in view of the apparent unification of all three gauge couplings in 
the MSSM in the absence of any intermediate scale. Similar mass patterns are, however, 
expected even if this GUT group first breaks to the Pati-Salam group. 

Our main goal is to compare and contrast the SUSY signals predicted by the nonuniversal 
scenario with those of the conventional one (D = 0) which has 10 roughly degenerate squarks 
of L and R types excluding the stop (t). We shall restrict our analysis to SUSY searches 
at various phases of the Tevatron. Special attention will be paid to the following issues: i) 
Is it possible to get a signal at Run II in some region of the parameter space (not already 
excluded by LEP searches) of the nonuniversal model with mo 3> mi/2, which corresponds 
to rriq ^> nig in the conventional scenario? ii) If signals are seen at the Tevatron, can they 
always be accommodated in both models, or might a distinction between the two be possible? 

In the next section we shall present the strategy of this analysis, the choice of the 
parameters and the observables. The results will be presented in section 3. Our conclusions 
and future outlook will be given in section 4. 



2 Strategy, the choice of parameters and observables 

In order to analyze signals for the production of squarks and gluinos in a systematic way the 
large nig — nig parameter space is usually divided into the following three broad regions: 
(A) nig ^> nig, (B) nig « m g and (C) nig <C m g . 

Out of these, region (C) is not accessible to first and second generation squarks in any 
grand desert scenario: even if their masses vanish at the GUT scale, large contributions 
from gluino-quark loops imply that at the weak scale nig ^ 0.85m g . However, some third 
generation squarks might be significantly lighter than the gluino. 

Our analysis will be based on the following strategy. We shall choose a set of GUT scale 
parameters viz., mo and m^, which evolves into weak scale squark and gluino masses in 
regions (A) and (B) in the conventional scenario. Moreover, mi/2 is chosen such that the 
resulting chargino mass is always consistent with the LEP bounds. As mentioned in the 
Introduction, we assume that the masses of b and f at scale Mq are the same as those of 
the corresponding first generation sfermions, i.e. we ignore the effect of RG scaling between 
A and Mq- Since we will concentrate on a moderate value of tan/?, this should be a good 
approximation for both scalar r states as well as for 6r. This assumption will probably 
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over-estimate the masses of third generation 577(2) doublet squarks slightly, but we expect 
this effect to be smaller than the new D— term contributions which are the main focus of 
our studyf] 

We have chosen tan/5 = 3 and the trilinear coupling A$ = 0. The latter choice of course 
has no direct bearing on the phenomenology of the first two generations, which is the focal 
point of our attention. Now for D = the magnitude of the Higgsino mass parameter 
fj, is fixed by the radiative electroweak symmetry breaking condition fllSl . We shall work 



throughout with negative /i, since otherwise LEP searches require m~ g to be too large for 
gluino searches at the Tevatron to be viable. 

In the presence of nonuniversal Higgs masses the magnitude of /i fixed from the EW 
symmetry breaking condition can in principle be quite different from that in the conventional 
scenario || . However, in order to compare our results with the predictions of the conventional 
scenario, we shall work with the same /i as in the D = case. This does not affect the cross 
sections for the production of the sparticles we are interested in, with the sole exception of 
the light scalar top eigenstate t\. Its mass depends on /i through the off-diagonal entries of 
the stop mass matrix. However, as will be discussed in more detail below, even a rather light 
t\ would not contribute very much to the signals we are interested in. We have therefore 
ignored its contribution completely; however, we have checked that the value of mi that 
follows from our (in this case not very reliable) assumptions satisfies the relevant search 
limits. 

The value of /x also affects the leptonic branching ratio of the second neutralino x\ ^ & 
decays into 3-body modes. The reason is that the Zx\x\ coupling only proceeds through the 
higgsino components of both neutralinos, and therefore scales roughly oc l//i 2 . As a result, 
slepton exchange contributions to X2 decay can compete with Z exchange contributions as 
long as slepton masses do not greatly exceed However, we will see that our model already 
allows considerable variation of this branching ratio even for fixed m\ii and m ; allowing 
H to vary away from the value predicted in the universal scenario would therefore not lead 
to new (combinations of) signatures, although it might shift the regions in parameter space 
where final states containing leptons are important. 

Having chosen the parameters, it is now straightforward to compute cross sections at 
the parton level in both scenarios. We shall consider the cross section o n corresponding to 
events with n leptons + jets + missing transverse energy (I$ T ) coming from all types of 
squark and gluino (excluding t) production channels. In particular o"o and 02 will turn out 
to be crucial observables in our analysis. We also include the cross section of hadronically 
quiet trilepton events 03; in our set of observables. These arise from xf — x® production and 
depend strongly on the slepton masses through the leptonic branching ratio of In this 
exploratory study we do not attempt any event simulation, nor do we apply the cuts required 
to extract these signals. We should clarify that we assume the experimental definition of the 
hadronic observables o n to include fairly hard cuts on the energies of the jets and/or on the 
missing transverse energy; this is why we do not include contributions from the production 
of electroweak gauginos here, even though they may give the biggest contribution to the 



*Note that we do not assume the Yukawa couplings to unify. In the simplest 50(10) model this assumption 
would require a large value of tan/3, leading to reduced masses for all third generation sfermions. Signals 
requiring the presence of energetic electrons or muons will then be greatly depleted even in the conventional 
scenario H. 
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total SUSY cross section at the Tevatron. The same hard cuts would also greatly reduce the 
contribution of a light t% to these final states. Moreover, we only count "primary" electrons 
and muons as "leptons"; we ignore leptons that may result from the decay of c— or b— 
quarks or r— leptons. 

In the absence of an event simulation which gives the efficiencies of the kinematical cuts 
to reduce the SM background, a full appreciation of the observability of the signals we have 
considered is not possible. Fortunately published results of such Monte Carlo simulations 
allow us to at least estimate the minimal cross sections (before cuts) that lead to detectable 
signals at future runs of the Tevatron. The present bounds, from data taken during Run I, 
can be gleaned from experimental searches. Both these searches and theoretical studies of 
the reach of Tevatron upgrades have generally been performed in the conventional scenario. 
However, we believe that the efficiencies should not change by more than a factor of 2 when 
allowing non- vanishing values of D, as long as we compare scenarios with similar values of 

Our estimates for the minimal detectable values of Oo, 02 and 03; are collected in Table 1. 
The limits on <7o and a^i from previous Tevatron runs have been extracted from experimental 
searches [[3], The future sensitivity in the jets plus missing E T channels (00,0-2) as well 
as in the tri-lepton channel (03;) has been estimated using results from refs.Q] and [17|, 
respectively. We also used ref.[|J to estimate the sensitivity of Run I to o 2 , since the only 
published experimental analysis of the jets plus di-lepton channel fl8j is based on just 19 
pb -1 of data.Q 



Signal 


Run I 


Run II 


TeV33 


o 


5 


0.9 


0.10 


o 2 


0.6 


0.14 


0.03 


032 


0.25 


0.03 


0.009 



Table 1: The minimal detectable values of the cross sections for jets plus missing Et final states with zero 
(cto) or two (02) electrons or muons, as well as for the relatively hadron-free three lepton final state (173;). 
The second, third and fourth column show the sensitivity from present data, from Run II (2 fb _1 at ^/s = 2 
TeV) and from TeV33 (25 fb^ 1 at ^fs — 2 TeV), respectively. All cross sections are in pb. 



3 The spectra and the signals 

The sfermion masses at the GUT scale can be parametrized by 0: 

m l L = m l R = m l R = m l + 0.5Dm 2 ; (la) 

OOOO / \ 

m dfl = m~ L = m — 1.5Dm , (lb) 

^Both CDF and D0 have announced ]l9| preliminary limits from searches for final states containing two 
leptons. However, the CDF search only covers the signal from like-sign di-leptons, which according to ref.Q 
has a somewhat smaller reach than the inclusive di-lepton sample, while the D0 search did not include 
events containing muons. 
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where the unknown parameter D can be of either sign. As stated earlier we use the same 
mo for the third generation, which amounts to neglecting the running of their masses above 
M G . 

Requiring positivity of the sfermion mass squares at the GUT scale imposes the restric- 
tion —2 < D < 0.66. Slightly stronger constraints on D arise from experimental searches. 
In particular, LEP searches ]20| for SU (2) doublet and singlet sleptons translate into upper 
bounds on \D\ for positive and negative D, respectively. For large mo, the LEP bound 
on the mass of the lighter sbottom eigenstate b\ gives the most stringent bound on \D\ for 
D < 0. 



D 





0.4 


0.6 


-0.75 


-1.25 




646.7 


702.4 


728.6 


526.6 


428.3 


m dr, 


650.7 


706.0 


732.2 


531.5 


434.2 




644.3 


700.2 


726.5 


523.7 


424.6 


m dn 


645.1 


436.7 


279.3 


915.8 


1058.5 


m H 


525.6 


592.7 


623.6 


367.9 


203.5 


m~ eL 


606.1 


376.7 


170.9 


888.8 


1035.2 


m- eR 


602.6 


662.1 


689.8 


471.5 


358.3 


m v 


601.8 


369.9 


155.2 


885.9 


1032.8 



Table 2: Squark and slepton masses in GeV at the weak scale for different values of D with mo = 600 GeV 
and mi/2 = 105 GeV. 

We first study the signals for the choice mo ~ 6my 2 , which corresponds to region (A) 
introduced in the last section. This is the for our purposes most interesting region, since 
here a signal at the Tevatron or its upgrades is not expected in the conventional scenario. 
Examples for sfermion masses in this region of parameter space are given in Table 2, for 
tuq = 600 GeV and mi/2 = 105 GeV and several values of D. The gaugino spectrum is 
practically independent of D: m- g = 332.1 GeV, M~± = 95.3 GeV, m~o = 95.5 GeV, and 
m~o = 46.5 GeV. 

Two points are worth noting : 

• The masses of SU(2) doublet sleptons II and cIr— type squarks can be reduced signifi- 
cantly for D > 0; for D > 0.55, g — > d R d 2-body decays become possible. 

• Due to the smaller coefficient of D in eq.(|la|) as compared to ( |I~b|) , the reduction of 
the masses of u^^ d^ and ur squarks as well as Ir sleptons for D < is at first less 
significant. However, for large negative D (D < —0.8) the bi becomes lighter than 
the gluino; in fact, as mentioned in the previous section, the lower bound on D comes 
from the experimental bound on m^ in this scenario.^ 

* A recent study [^2| of the sparticle spectrum predicted by the minimal 5O(10) SUSY GUT also empha- 
sized the reduction of mg relative to the other squark masses. However, as mentioned in Sec. 2, minimal 
5O(10) requires large values of tan/3, leading to a large bottom Yukawa coupling. Successful electroweak 
symmetry breaking then requires D > 0. As a result, the light bi found in ref. |^2| is mostly an SU(2) singlet, 
while our more general scenario allows D < and a light b^. 
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These changes of the scalar spectrum lead to intriguing features of gaugino branching 
ratios. In the electroweak gaugino sector, BR(%2 —* Xi^) increases significantly for D > 0, 
due to the reduction of mj . For example, we have BR(%2 — * Xi^) — 0.02, 0.07 and 0.17 
(for each lepton flavor) for D = 0, 0.4 and 0.6, respectively. As shown in Fig. 1, for D ;$ 0.5 
this results in an appreciable enhancement of the dilepton signal (cr 2 ) compared to the 
conventional scenario. The clean trilepton signal (173;) grows monotonically with increasing 
D. 

Reducing D below zero at first leads to little change of BR(x2 — > Xi^)- Note that in 
the conventional scenario an increase of m beyond 600 GeV results in an increase of this 
branching ratio, due to reduced destructive interference between Z and slepton exchange 
contributions. In our case reducing D below zero leads to lighter SU (2) singlet sleptons, but 
heavier doublet sleptons; the two effects appear to largely cancel in the leptonic branching 
ratio of x%- However, for D < — 1.2 the reduced b L mass leads to a rapid increase of 
BR(x2 ~* Xipb), and a corresponding reduction of the leptonic branching ratio, as well as of 

&31- 

The chargino decays do not exhibit any remarkable feature over the entire region of D 
considered above, since they are dominated by virtual W exchange contributions. 

The gluino branching ratios for D = and 0.4 are:[] 
BR(g -> XiQQ') = 0-47, 0.37 (qq' = ud, cs), 
BR{g-*x\qq) = V.l^ 0.35 (q = u, d, c, s, b), 
BR(g -> xlqq) = 0.39, 0.28 (q = u, d, c, s, b). 

Note that the branching ratio for g — > XiQQ more than doubles when going from D = to 
D = 0.4. This is a consequence of the reduced dn masses; recall that SU(2) singlet squarks 
couple only very weakly to X2 an d Xi-, which are predominantly 577(2) gauginos. This leads 
to an increase of do (see Fig. 1), a hardening of the missing E? spectrum, and an increase 
of the fraction of gluino pair events containing b— jets. 

These changes become even more dramatic for D > 0.55. Now d^— type squarks become 
lighter than the gluino, so that all gluinos decay to d^s with BR = 1. Almost all the d^s in 
turn decay directly into the LSP. This suppresses 02 severely compared to the conventional 
case. In addition, gdn production as well as the pair production of dR— type squarks lead to 
a further increase of a , by up to a factor of 2. Recall that the clean trilepton signal (<t 3 ;) is 
also enhanced relative to the universal case. This combination of significant o"o and 173; but 
nearly vanishing 02 cannot be accommodated in the conventional scenario for any choice of 
the parameters. 

As mentioned above, taking D < decreases the masses of SU(2) doublet as well as 
ur— type squarks. This leads to increased destructive interference between s— and t— channel 
contributions to g pair production; as a result, the corresponding cross section drops by about 
a factor 1.6 when D is decreased from to —0.9. We also saw that the decrease in squark 
mass is especially pronounced for 6^, whose mass is also reduced by RG-running from the 
GUT to the weak scale. Reducing D below therefore also increases the fraction of g 
decays that produce b— jets; we saw a similar effect for D > due to the reduction of vn~ hR . 

^We ignore a small contribution from g — > X^tb. For our choice of gluino mass the corresponding partial 
width is suppressed by phase space. Moreover, its exact value depends sensitively on the masses and mixing 
angle of the stop squarks, which in turn depend on details of the GUT model, as discussed earlier. 
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However, while this fraction saturates at 1/3 for D > 0, where all three cIr— type squarks 
become light, it reaches 100% for D < —0.9, where g — ■> b^b is the gluino's only 2-body 
decay mode. Nearly all bi in turn decay into xV 3 - Since leptonic \2 decays always lead to 
lepton pairs with opposite charges, the fraction of like-sign di-lepton pairs in <ji becomes 
very small here. 

For D < - 1, b L pair production begins to contribute significantly to o~o and cr 2 - For 
example, for D = —1.3 (—1.4) we have = 174 (115) GeV, giving a(b L b* L ) ~ 1.0 (8.5) 
pb. Due to the presence of (at least) two hard b— jets in the event, a dedicated sbottom 
pair search is expected to be sensitive to much smaller cross sections than more generic 
SUSY searches. A recent analysis concludes that data from from Run II (TeV33) should 
allow to search for b^ with mass up to at least 185 (205) GeV even if all X2 decay into 
light quarks; in the case at hand this would cover the range D ^ — 1.25. Scenarios with 
D ^$ — 1.4 can perhaps even be probed with existing data from Run I.[] Recall also that 
for large negative D, BR(x2 ~> Xi^b) is enhanced, leading to events with up to 6 6— jets! 
In the conventional scenario one would expect a and a 2 to be dominated by sbottom pair 
production only at very large values of tan/3; and even there the b\ — g mass splitting is less 
than for large negative D in the non-universal case. Moreover, in the conventional scenario 



one also expects final states to be rich in r— leptons ||14| : this need not be the case in the 
non-universal scenario. 

It is clear from Table 1 and Figure 1 that if rriq 3> m g , the latter being tightly constrained 
by LEP bound on M~±, the conventional scenario offers very little hope for a discovery in 
any of the three channels up to the TeV33 phase. In contrast, the non-universal scenario 
allows for cross sections in the striking range of Run II experiments, from sbottom pair 
production for large negative D, and in the tri-lepton channel at large positive D. Unfor- 
tunately observation of a signal in the tri-lepton channel by itself would not be sufficient 
to discriminate between the two scenarios. In the case at hand this distinction should be 
possible using TeV33 data, since there cxo and a^i become measurable over the entire range 
of D shown in Fig. 1. As noted earlier, the combination of large <To and a^i with small 02 
would exclude the conventional scenario. Additional clues can be obtained by correlating 
kinematical distributions, which are sensitive to sparticle masses, with measurements of the 
total cross sections. 

The nonuniversal model allows for an observable trilepton signal even for higher m gi 
albeit for much lighter sleptons (D ^ 0.6). For example, with D = 0.6, m 1( / 2 = 130 GeV {m g 
= 420 GeV), mj L = 190 GeV, a 3 i = 0.09 pb should be easily detectable already at Run II. 
The corresponding cr is, however, unobservable even at TeV33. 

We next consider region (B) of the parameter space. We shall work with m 1( / 2 = 105 
GeV and mo = 300 or 200 GeV. This is because even the phenomenology of the conventional 
scenario depends sensitively on mo. Additional comments will be made on higher m g as and 
when appropriate. For itlq = 300 GeV, the following mass patterns are obtained: 



*A search by the D0 collaboration @ was able to exclude m gi < 115 GeV if BR(bi -> x?6) ~ 100% 
and m~ ^ 20 GeV; neither of these conditions is satisfied in our case. However, during Run 1 CDF had 
significantly better b— tagging capabilities than the D0 detector. 



D 





0.5 


-1.00 


-1.75 


rriu L 


392.6 


421.5 


327.4 


268.3 


m,- 

d R 


389.9 


285.6 


541.2 


631.3 


m ~br. 


339.9 


372.3 


261.1 


181.6 




312.2 


164.3 


488.2 


586.5 


VTiy 


303.9 


147.9 


482.9 


582.1 




305.3 


341.6 


215.0 


105.0 



Table 3: Squark and slepton masses (in GeV) at the weak scale for different values of D with too = 300 GeV 
and mi/2 = 105 GeV. 



The spectrum for the gaugino sector is: m g = 320 GeV (the difference from the previous 
case arises due to radiative corrections to the gluino pole mass from squark-quark loops), 
M~± = 95.16 GeV, m~o = 95.06 GeV, and m~ = 46.71 GeV. 
Points to be noted are : 

• In the conventional scenario the squarks are somewhat heavier than the g. 

• Increasing D from zero again reduces and m^; in particular, for D > 0.5, mj fl < 
trig for all three generation of dn— type squarks; and 

• For large negative D (D < — 1) the hi again becomes lighter than the gluino, but the 
theoretical lower bound on its mass is higher than in the previous case, because of the 
LEP bounds on m^. 

In contrast to the previous case, all possible squark and gluino production channels can 
contribute to a and cr 2 - As shown in Fig. 2, the size of these cross sections will be typically 
larger than the conventional case even for moderate magnitudes of D irrespective of its sign. 
o"o still stays below the current bound of about 5 pb over the entire range of D; however, the 
region with D ^ — 1 should be accessible in this channel at Run II, and the entire range of 
D can be probed at TeV33. 

Unlike in the scenario with large mo we now observe a significant increase of the leptonic 
branching ratio of x.2 f° r large negative D. Note that now slepton exchange contributions 
dominate the leptonic partial width of even in the conventional scenario. Decreasing D 
away from zero therefore at first actually reduces the leptonic branching ratio of x\i due to 
the increase of m^; note that for roughly equal mass the exchange of SU{2) doublet sleptons 
contributes more, due to their stronger coupling to x%- However, for D ^$ — 1 the reduction 
of mg R over-compensates the increase of mg L , leading to an increase of 03/; in fact, the region 



around D = —1.75 might already be excluded by existing searches [1C] for tri-lepton events. 
For even more negative values of D, the decay Xi ~ * ^1 + v t opens up. Since we do not 
include secondary leptons from r decays when computing 03/, this cross section essentially 
vanishes once that decay mode dominates. In reality a measurable tri-lepton signal may 
survive, if softer cuts on the leptons are employed [[171] ; in addition, one may look for events 
with two leptons and a r— jet |14] . Note that unlike xt, X2 a ^ so couples to and /i^; the 
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leptonic branching ratio of X2 therefore remains large all the way out to the LEP-imposed 
lower bound on D. 

Together with the reduction of mg this can lead to promising signals from events with 
at least one hard b— jet and two opposite-charge, same-flavor leptons; this might extend the 
reach of Run II to smaller values of |Z)| than those accessible to sbottom searches based 
on final states without leptons p3| . Even these kinds of searches should be able to probe 
D < — 1.75 (—1.6) at Run II (TeV33). Moreover, even a moderately negative D leads to 
a significant increase of Oi- Recall that in this region all SU(2) doublet squarks, as well 
as Mr— type squarks, are somewhat lighter than in the conventional scenario (for the same 
value of mi/2 and mo), leading to a significant increase of the total cross section for the 
production of strongly interacting sparticles. This over-compensates the modest decrease 
of the leptonic branching ratio of for moderately negative D. Nevertheless this region of 
parameter space can be probed in the inclusive di-lepton channel only at TeV33. 

Increasing D from zero leads to quickly increasing leptonic branching ratio of x%i an d a 
corresponding rise of a^; as for negative D, the region just below the largest value of \D\ 
allowed by slepton searches at LEP is again at best marginally consistent with current search 
limits for tri-lepton final states. For the largest LEP-allowed value of D, the decays x° ~* vv 
open up. The sneutrino then decays into the invisible channel v — > y§y (the so-called VLSP 
scenario [p5fl), leading to a vanishing u 3; . Apart from the most extreme regions, the entire 
range of D can easily be covered in the tri-lepton mode at TeV33, and much of it should 
already be accessible to Run II experiments. Finally, sizably positive D again suppresses 02, 
first due to the dominance of g — > d^d decays (for D > 0.35), and, for the largest allowed 
value of D, due to X2 ~* vv decays. 

For itlq = 200 we obtain the following sfermion spectra: 



D 





0.4 


0.6 


-0.75 


-1.75 


m UL 


326.6 


339.1 


345.2 


301.7 


264.4 


mi 

an 


323.4 


282.0 


258.8 


389.3 


462.7 


m h 


293.6 


307.5 


314.2 


265.6 


222.9 


m~ eL 


218.0 


149.9 


99.7 


307.4 


396.3 


me R 


207.8 


227.0 


236.1 


166.0 


81.7 


mc, 


205.9 


131.7 


69.5 


298.9 


389.8 



Table 4: Squark and slepton masses (in GeV) at the weak scale for different values of D, for too = 200 GeV 
and TO1/2 = 105 GeV. 

The gaugino masses are nearly the same as for m = 300 GeV. The behavior of our three 
cross sections, shown in Fig. 3, is qualitatively similar to the previous case. However, due 
to the reduced importance of the D— terms as compared to the universal contribution from 
gluino-quark loops, even in the non-universal scenario the mass splitting between different 
first or second generation squarks now amounts to at most 50%. As a consequence the 
changes of the three cross sections o"o, 02 and cr 3 ; with varying D are somewhat smaller than 
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for larger values of m . In particular, a now increases by less than a factor of 3 as D is 
decreased from zero to its most negative allowed value. On the other hand, owing to the 
reduction of the overall squark mass scale, the entire range of D can now be probed in the 
jets plus missing Et channel at Run II of the Tevatron. 

Somewhat larger variations can occur for the leptonic cross sections. Since sleptons 
are still significantly heavier than \\ an d Xi in the conventional scenario (D = 0), the 
leptonic branching ratios can change considerably when D is allowed to take values close to 
the present upper or lower bounds. Due to the reduction of the slepton masses compared 
to the previous cases, these cross sections are also generally larger than for scenarios with 
larger m . In particular, the tri-lepton signal should be detectable over the entire range of 
D at Run II. The availability of at least two distinct SUSY signals should make it easier 
to distinguish between the universal and non-universal scenarios; on the other hand, the 
expected differences in the squark spectrum are smaller than before, as already noted. 

Yet smaller values of mo continue the trend of going from 300 to 200 GeV. The maximal 
possible mass splitting between squarks keeps decreasing even in the non-universal scenario. 
Significant differences in the slepton spectrum remain possible; as long as sleptons are heavier 
than SU (2) gauginos in the conventional scenario, these changes of slepton masses affect the 
leptonic branching ratios of electroweak gauginos. However, for m < 85 GeV, these gauginos 
can decay into on-shell sleptons even in the universal scenario; their leptonic branching ratios 
are then much less sensitive to D. Moreover, LEP searches severely limit the maximal allowed 
value of D for such small values of m . 

4 Summary and conclusions 

In this paper we have studied the impact of 5*0(10) D— term contributions to scalar masses 
on the search for supersymmetric particles at the Tevatron collider. These contributions have 
opposite signs for sfermions residing in the 5* representation of SU(5) and those residing in 
the 10 representation. Any non-zero value of D will therefore reduce the masses of some 
sfermions, while making others heavier. This usually results in an increase of the total cross 
section for sparticle production at hadron colliders. Heavy sparticles simply decouple, so 
that further increases of their masses have little impact on any observables; however, the 
cross section for the production of the lighter sfermions keeps increasing as their masses are 
reduced. The only exception to this rule occurs for large m and moderately negative D, 
where the reduction of most squark masses actually reduces the gluino pair cross section 
through interference effects. However, we found that this can reduce the total cross section 
for the production of strongly interacting sparticles by no more than a factor of 2; in contrast, 
for extreme (but allowed!) values of D this cross section can be enhanced by more than an 
order of magnitude compared to models with degenerate squarks. 

Similarly, non-degeneracy in the slepton sector tends to increase leptonic branching ratios 
of electroweak gauginos, in particular those of leading to an increase of the relatively 
hadron-free tri-lepton cross section from xtxf production. For large negative D this cross 
section can be sufficiently large to yield unambiguous signals at the next run of the Tevatron 
collider even if the average sfermion mass lies well above the gluino mass, which in turn lies 
well above its lower limit derived from chargino searches at LEP (assuming gaugino mass 
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unification). However, if D is very close to its LEP-imposed upper or lower limit, the tri— 
lepton signal is greatly depleted, since there decays X2 ~^ Xi^ u an d Xi ~^ TxV T1 respectively, 
become possible. 

One interesting feature of the non-universal scenario is that one can find combinations 
of parameters where the tri-lepton and multi-jet plus missing Et signals are both large, but 
the jets plus di-lepton plus missing Et signal is very small. This last signal results from 
cascade decays of the gluino and SU(2) doublet squarks, which are suppressed for D > 0, 
where cLr— type squarks can be the only light squarks. Such a combination of signals cannot 
occur if all sfermions are degenerate at the GUT scale. 

Models with non-vanishing D also introduce new signals, which do not hold much 
promise if squarks are approximately degenerate. These signals originate from the copi- 
ous production of light sbottom squarks. In models with small or moderate tan/? studied 
here, light b^ squarks can occur for negative D. Since b\ dominantly decays into x^b, b\ pair 
production (as well as gluino production followed by g — > b\b decays) gives rise to several 
signals. Since for large negative D the branching ratio for %\ — > x\bb decays can be en- 
hanced, searches for final states containing 4, or even 6, b— jets may offer the best reach. For 
somewhat smaller \D\ or smaller mo, searches for final states containing two hard b— jets and 
two opposite-charged, same-flavor leptons might hold even more promise. This reinforces 
the conclusions of ref.|22| regarding the importance of SUSY signals containing hard b— jets. 

The maximal possible mass splitting in the squark sector becomes smaller for smaller 
ratios of the GUT scale parameters mo/m^, due to the increased importance of flavor- 
universal contributions to squark masses from gluino-quark loops. For m < 2mi/ 2 the main 
impact of sfermion non-universality on SUSY signals at the Tevatron comes from changes of 
the leptonic branching ratios of electroweak gauginos, which remain sensitive to D as long as 
mo^mi/2. For even smaller values of mo/m^j evidence for non-universality in the sfermion 
sector can probably best be found in studies of slepton pair production. However, these may 
have to be postponed until a lepton supercollider has been built. 

In summary, in 5*0(10) models with sizable D— term contributions to scalar masses, 
SUSY searches at future collider runs of the Tevatron hold even more promise than in the 
conventional mSUGRA model. New regions of the (mo,m 1; /2) plane become accessible to 
"standard" searches involving jets and/or leptons plus missing transverse momentum, and 
searches for new final states containing hard b— jets may increase the reach even further. 
Under favorable circumstances the pattern of observed signals will allow to discriminate 
between universal and non-universal scenarios. This once again underscores the importance 
of searching for supersymmetry in as many different channels as possible. 
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Figure 1: Cross sections (in pb, without any cut) for different final states at the Tevatron 
collider with y/s = 2 TeV as a function of D. The other relevant parameters are = 105 
GeV, m = 600 GeV and tan/5 = 3, while /i is fixed by imposing the radiative electroweak 
symmetry breaking condition for D — 0. 




Figure 2: As in Fig. 1, but with m = 300 GeV. 
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Figure 3: As in Fig. 1, but with m = 200 GeV. 
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